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A search for depletion of the combined flux of active neutrino species over a 735 km baseline 
is reported using neutral-current interaction data recorded by the MINOS detectors in the NuMI 
neutrino beam. Such a depletion is not expected according to conventional interpretations of neu- 
trino oscillation data involving the three known neutrino flavors. A depletion would be a signature 
of oscillations or decay to postulated noninteracting sterile neutrinos, scenarios not ruled out by 
existing data. From an exposure of 3.18 x 10 20 protons on target in which neutrinos of energies 
between ~500 MeV and 120 GeV are produced predominantly as i/ M , the visible energy spectrum of 
candidate neutral-current reactions in the MINOS far-detector is reconstructed. Comparison of this 
spectrum to that inferred from a similarly selected near-detector sample shows that of the portion 
of the Vp, flux observed to disappear in charged-current interaction data, the fraction that could be 
converting to a sterile state is less than 52% at 90% confidence level (C.L.). The hypothesis that 
active neutrinos mix with a single sterile neutrino via oscillations is tested by fitting the data to 
various models. In the particular four-neutrino models considered, the mixing angles #24 and #34 
are constrained to be less than 11° and 56° at 90% C.L., respectively. The possibility that active 
neutrinos may decay to sterile neutrinos is also investigated. Pure neutrino decay without oscilla- 
tions is ruled out at 5.4 standard deviations. For the scenario in which active neutrinos decay into 
sterile states concurrently with neutrino oscillations, a lower limit is established for the neutrino 
decay lifetime r 3 /m 3 > 2.1 x 10~ 12 s/eV at 90% C.L. 

PACS numbers: 14.60. St, 12.15.Mm, 14.60.Pq, 14.60. Lm, 29.27.-a, 29.30.-h 



I. INTRODUCTION 



Compelling evidence has been presented demonstrat- 
ing the disappearance of muon and electron neutrinos 
as they propagate from their production point. Dis- 
appearance of muon neutrinos has been observed from 
neutrino fluxes originating in the atmosphere [l], Q and 
accelerator beams [!, 0| ; disappearance of electron neu- 
trinos has been observed with neutrino fluxes from the 
Sun 0, Q and terrestrial reactors Q . Super-Kamiokande 
and other atmospheric neutrino experiments were the 
first to report significant deficits of charged-current 
interactions from neutrinos propagating over baselines 
larger than several hundred kilometers. The K2K and 
MINOS experiments have observed the disappearance 
using accelerator-beam neutrinos propagating over fixed 
baselines of 250 km and 735 km, respectively. All experi- 
ments reporting muon-neutrino disappearance favor pure 
Vp — > v T oscillations as the explanation for the observed 
disappearance of [3I. fsj-TTol] . However, more exotic sce- 
narios in which one or more sterile neutrinos, v a , mix 
with the three active neutrino species remain as viable 
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alternatives. 

Long-baseline experiments such as MINOS provide 
an opportunity to test alternative scenarios by com- 
paring the observed neutral-current interaction rates 
in near and far-detectors. Since all active neutrinos, 
v e , i/p, and v T , participate in the neutral-current inter- 
action, this comparison provides a sensitive probe to the 
existence of processes that deplete the flux of active neu- 
trinos between the two detectors. If neutrinos only oscil- 
late among the active flavors, the rate of neutral-current 
interactions at the far site of a long-baseline experiment 
remains unchanged from the null-oscillation prediction. 
However, if another process occurs concurrently with ac- 
tive neutrino oscillations, the rate of neutral-current in- 
teractions at the far site may be different. Two such pos- 
sibilities have attracted considerable attention and are 
the focus of the analysis reported here, (i) active neutri- 
nos oscillating to v Sl and (ii) neutrino decay in conjunc- 
tion with oscillations. 

The possible existence of one or more sterile neutrinos 
that do not couple to the electroweak current but do mix 
with the active flavors has been discussed extensively in 
the literature 

GH3. 

The existence of sterile neutrinos 
would provide new degrees of freedom that could help 
clarify certain outstanding problems with the neutrino- 
mass spectrum [hi ] and with heavy element nucleosyn- 
thesis in supernovae (l5| . A recent search for neutrino 
oscillations in a short baseline experiment provides no 
evidence for transitions that would imply existence of 
sterile neutrinos flfjj . 

The coupling between sterile neutrinos and the ac- 
tive neutrinos would likely involve the third mass eigen- 
state. Observations by the SNO experiment show the 
total flux of active neutrinos from the Sun to agree with 
expectations from solar models @, thereby limiting the 
potential coupling of the first or second neutrino-mass 
eigenstates to a sterile neutrino. Additionally, the Super- 
Kamiokande experiment strongly disfavors pure — > v s 
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mixing [8J, |9[ , but does not exclude an admixture of sub- 
dominant i/„ — > v s mixing with the dominant v„ — > v T 
mixing. MINOS has recently carried out the first ded- 
icated search at fixed long-baseline for oscillating to 
both v T and v s (l7| . The analysis presented here uses a 
larger exposure and extends the earlier analysis by con- 
sidering specific models in which a sterile neutrino state 
is incorporated into the neutrino mixing matrix. 

The possibility that a neutrino may decay into a sterile 



state |18[ is also explored in this work. In this scenario, 
the mass eigenstate 1/3 is unstable and allows active neu- 
trinos to decay into undetectable final states. The decays 
would give rise to an anomalous depletion of the neutral- 
current event rate observed at the far-detector. The oc- 
currence of pure neutrino decay, without oscillation, has 
already been shown by MINOS and Super-Kamiokandc 
to be highly disfavored [T3, Ell ■ The analysis reported 
here represents the first direct test of the neutrino-decay- 
with-oscillations scenario in a long-bascline experiment. 



II. NUMI BEAM AND MINOS DETECTORS 



The MINOS near-detector is located 1.04 km 
downstream of the target, has a total mass of 
980 metric tonnes, and lies 103 m underground at 
Fcrmilab. The detector consists of two sections, a 
calorimeter encompassing the upstream 121 planes and a 
spectrometer containing the downstream 161 planes. In 
both sections, one out of every five planes is fully covered 
with 96 scintillator strips attached to the steel plates. In 
the calorimeter section, the other four out of five planes 
are partially covered with 64 scintillator strips, whereas 
in the spectrometer section no scintillator is attached to 
the steel. The far-detector is 734 km downstream of the 
near-detector, has a total mass of 5400 metric tonnes, 
and is located in the Soudan Underground Laboratory, 
705 m below the surface. It is composed of 484 fully 
instrumented planes organized in two supermodules |4|. 
The fiducial masses used for the near and far-detectors 
are 27 metric tonnes and 3800 metric tonnes respectively. 
The near-detector steel is magnetized with an average 
field intensity of 1.3 T whereas the far-detector has an 
average field of 1.4 T in the steel. 



Neutrinos from the NuMI (Neutrinos from the Main 
Injector) beam [20| originate from decays of pions and 
kaons produced in the beamline target; a significantly 
smaller contribution arises from subsequent muon decays. 
The secondary mesons are created using 120 GeV pro- 
tons extracted from the Fermilab Main Injector incident 
on a graphite target. The proton extraction occurs in 
10 /is spills with a 2.2 s cycle. Positioned downstream of 
the target are two parabolic magnetic horns which fo- 
cus 7r + and K + secondary particles. The focused mesons 
proceed into a 675 m long evacuated decay pipe, where 
they may decay into muons and neutrinos. The remnant 
hadrons are stopped by a beam absorber placed at the 
end of the decay pipe. The tertiary muons are stopped 
by 240 m of rock between the end of the decay volume 
and the near-detector cavern so that only neutrinos reach 
the near-detector. The neutrino energy spectrum can be 
changed by adjusting cither the horn current or the po- 
sition of the target relative to the horns. The data em- 
ployed in this analysis were obtained using the low-energy 
beam configuration, in which the peak neutrino energy 
is 3.3 GeV [J, and correspond to a far-detector exposure 
of 3.18 x 10 20 protons on target, collected during the pe- 
riod of May 2005 to July 2007. In this configuration, ac- 
cording to Monte Carlo simulations, the neutrino flavor 
composition of the beam is 91.8% v^, 6.9% F^, and 1.3% 
v e + V e . For this analysis the neutrinos and antineutrinos 
are assumed to oscillate with the same parameters. 

The MINOS detectors are planar steel/scintillator 
tracking calorimeters (2(j. The vertically oriented de- 
tector planes are composed of 2.54 cm thick steel and 
1 cm thick plastic scintillator. A scintillator layer is com- 
posed of 4.1cm wide strips. Each strip is coupled via 
a wavelength-shifting fiber to one pixel of a multianode 
photomultiplier tube (PMT) [2l|,[22|. 



III. DATA SELECTION 

A neutrino interacting in one of the MINOS detectors 
produces either a charged-current event with a charged 
lepton plus hadrons emerging from the event vertex or a 
neutral-current event with hadrons but no charged lep- 
ton in the final state. In either case, the particles in 
the final state deposit energy in the scintillator strips, 
which is converted to light and collected by optical fibers 
and converted to electronic signals by PMTs. The MI- 
NOS reconstruction algorithms use event topology and 
the recorded time stamps of the strips where energy was 
deposited to identify neutrino events inside the detector. 
Events must have at least four strips with signal to be 
considered in the analysis. Individual scintillator strips 
are grouped into either reconstructed tracks or showers, 
and the tracks and showers are combined into events Q ■ 
The vertex of each event is required to be sufficiently far 
from any edge of the detector to ensure that the final- 
state hadronic showers are contained within the instru- 
mented regions of the detectors. On average, each GeV 
of energy deposition in a neutral-current event induces 
activity in 12 scintillator strips. 

The Monte Carlo simulation of the neutrino beam uti- 
lizes FLUKA05 [23[ to model hadroproduction in the 
NuMI target and a GEANT3 [H| simulation of the NuMI 
beam line to propagate the particles exiting the target. 
The neutrino interactions in the MINOS detectors are 
modeled by the NEUGEN-v3 [H] program. The simu- 
lated neutrino flux is constrained to agree with the neu- 
trino energy spectra measured in the near-detector for 
nine different beam configurations Q. This procedure 
reduces the uncertainties due to the neutrino flux in the 
far-detector prediction. 
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A. Data integrity 

All of the data accepted for this analysis must pass a 
series of requirements on the beam and detector perfor- 
mance. The beam is required to strike within 2 mm of 
the center of the upstream face of the NuMI target, a 
segmented rectangular graphite rod measuring 6.4 mm 
in width, 15 mm in height and 940 mm in length 
The full width at half-maximum of the beam at the 
target is required to be between 0.1mm and 2.0 mm in 
the transverse horizontal direction and between 0.1mm 
and 1.5 mm in the transverse vertical direction. The 
minimum allowed beam intensity during a beam spill is 
0.5 x 10 12 protons on target. 

For all data used in this analysis, all detector subsys- 
tems that affect data quality are required to be in normal, 
stable modes of operation. Checks are made to ensure 
the coil currents that energize the magnetic fields are 
at their nominal values in both detectors. The timing 
between the detectors is synchronized using the Global 
Positioning System to within 1 fis to ensure that correct 
beam extraction timing is provided to the far-detector 
electronics. The timing synchronization is not affected 
by the timing resolution of the detectors, which is 18.8 ns 
and 1.9 ns for the near and far-detectors, respectively. 
The high voltage supplied to each PMT is required to be 
at its nominal value. 



B. Fiducial requirements 

Only those events reconstructed in the fiducial volume 
are included in the analysis. In both detectors the re- 
constructed event vertex is required to be more than 
50 cm from the edges of the instrumented regions and, 
in the case of the far-detector, more than 45 cm from the 
center of the magnetic coil that runs through the mid- 
dle of each detector plane. In addition, a longitudinal 
veto region comprising cither 30 planes at the front of 
the near-detector or four planes at the front of each of 
the far-detector supermodules eliminates events that en- 
ter through the first plane of a detector but may have 
originated outside the detector volume. To ensure good 
shower containment, the event vertex is required to be 
reconstructed more than 1 m from the last plane of each 
far-dctcctor supcrmodule and more than 1 m from the 
last plane in the calorimeter region of the near-detector. 
The sparsely instrumented spectrometer region of the 
near-detector is not included in the fiducial volume. The 
primary vertex of an event is assigned according to the 
event's reconstructed track vertex. However, for events 
without a reconstructed track, the vertex of the hadronic 
shower is used as the event vertex. 



C. Near-detector event selection 

The reconstruction algorithms are designed to han- 
dle the high rate of interactions occurring in the near- 
detector during running with typical intensities of 2.2 x 
10 13 protons on target per beam spill. At this intensity, 
an average of 16 neutrino interactions occur in the near- 
detector for each spill. For the majority of events the 
algorithms perform very well. However, for certain event 
subclasscss, shortfalls have been identified and quantified 
using special studies including low intensity beam data 
and visual scanning. Monte Carlo studies show that 8% 
of all reconstructed events classified as neutral-current 
interactions are assigned a value of recontructed energy 
that is less than half of the actual deposited energy in 
the detector for the simulated interaction; the remaining 
energy has been reconstructed as a separate event, re- 
sulting in an overestimate in the number of reconstructed 
events. In particular, the number of events with visible 
energy below 1 GeV that are classified as neutral-current 
candidates is overestimated by 34%. Therefore, these 
poorly reconstructed events affect the neutral-current en- 
ergy spectrum and care has been taken in identifying and 
removing them from the analysis. 

Reconstruction failures are classified into three main 
categories: (i) split events, (ii) leakage events, and (iii) 
incomplete events. Split events occur when a single neu- 
trino interaction results in two or more reconstructed 
events. Leakage events are due to incorrectly assigned 
event vertices causing neutrino interactions outside the 
fiducial volume to be reconstructed within it. The incom- 
plete event category is a looser classification that refers 
to further types of failures in shower reconstruction to 
be described below. In all three categories, the visible 
energy of a neutrino candidate may be underestimated, 
resulting in a background to neutral-current events at low 
energies. As near-detector data are used to predict the 
expected spectra at the far-detector, reconstruction fail- 
ures specific to the near-detector must be minimized. A 
set of selection requirements was developed [26[ to reduce 
the occurrence of these failure modes. After applying the 
requirements, detailed below, simulations show that the 
background of poorly reconstructed events having visible 
energy below 1 GeV has been reduced to 8%. 

Split events lead to double-counting of neutrino in- 
teractions and incorrect energy reconstruction. If two 
reconstructed events are caused by the same neutrino 
interaction, they are expected to appear close in time 
and in space. The time for each event has been taken 
to be the median of the recorded signal arrival times 
for event strips contained within five planes of the event 
vertex. To minimize the occurrence of split events, the 
time separation between an event and the closest other 
event in time is required to be |Ai| > 40 ns, as shown in 
Fig. [TJ A requirement that the spatial separation be- 
tween events along the beam direction |Az| > lm if 
40 ns < | At | < 120 ns is also employed to further elimi- 
nate split events. 
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FIG. 1: The distribution of time separation between events At 
for data (solid points) and Monte Carlo simulated data (solid 
histogram). A background component arising from poorly re- 
constructed events (hatched histogram) is confined to a nar- 
row region of low At values. Events accepted for further anal- 
ysis are indicated by the arrows. 



FIG. 2: Section of the distribution of the number of strips 
with nonzero pulse-height, per event, after all other selec- 
tions are applied. The region displayed, in which the contri- 
bution from poorly reconstructed events (hatched histogram) 
is significant, corresponds to low strip counts and represents a 
small fraction of the total number of events. The event range 
accepted for the analysis is identified by the arrow. 



Leakage events are typically cosmic-ray muons causing 
steep showers with a high concentration of hit strips in 
a small number of detector planes. These events can be 
removed by selecting on this topological characteristic. 
Thus, a requirement is placed on the ratio of the average 
number of active strips per plane to the total number of 
planes with active strips in the event, represented more 
concisely as (active event strips)/(active event planes) 2 . 
Only events for which the ratio is less than 1.0 are ac- 
cepted by the analysis. 

Another type of leakage event is due to secondary par- 
ticles from interactions occurring outside of the fiducial 
volume. In the partially covered planes of the near- 
detector, the steel is instrumented with scintillator to 
within 16 cm from the left-hand-side edges of the steel 
plate and 1.4 m from the right-hand-side edges as viewed 
from along the beam direction. Consequently, secondary 
particles may enter the detector laterally due to the 
sparse instrumentation on the sides. For such cases the 
reconstruction algorithm is likely to fail in associating 
hits to events. Nevertheless, the extra activity at the 
edges of the fully covered planes is recorded and can be 
used to veto events within a time window. The veto crite- 
rion uses the number of active strips and the pulse-height 
in the edge regions of the detector. An event is accepted 
if the number of active strips in the veto region recorded 
within a ±40 ns window around the event vertex time is 
less than four and the pulse-height deposited in the veto 



region is less than 2 MIP a . These veto criteria are applied 
to events with visible energy less than 5 GeV in which the 
number of planes assigned to the reconstructed shower is 
greater than the number of planes assigned to the re- 
constructed track, as leakage events are reconstructed as 
low-energy showers without a clearly defined track. 

Incomplete events arise when the shower reconstruc- 
tion fails to assign all event-related strips to the shower. 
This type of reconstruction failure occurs if there are 
large gaps in a shower or if the shower is generally sparse. 
In a majority of these cases, events have a very low num- 
ber of reconstructed strips. Figure [5] shows the distribu- 
tion of the number of strips for the events, after applying 
the selection requirements. To minimize the number of 
incomplete events in the near-detector data sample, an 
event is required to have total number of strips greater 
than four. 

In summary, the selection criteria applied to the near- 
detector data are as follows: (i) the modulus of the time 
separation between events, |Ai|, must exceed 40ns; (ii) 
if 40 ns < |At| < 120 ns, the modulus of the spatial 
separation between events, \Az\, must exceed 1 m; (iii) 
the ratio (active event strips)/ (active event planes) 2 must 
be less than unity; (iv) for events with less than 5 GeV 
of reconstructed energy in which the number of planes 
is larger in the reconstructed shower than in the recon- 



a Minimum ionizing particle, equivalent to the response produced 
by a 1 GeV muon traversing a detector plane at normal incidence. 



6 



structed track, the number of event strips reconstructed 
in the detector's veto regions should be less than four and 
the total pulse-height in those regions must be less than 
2MIP; and (v) the total number of strips reconstructed 
in the event must be more than four. Only events that 
satisfy these criteria are used for further analysis. 



D. Far-detector event selection 

In contrast to the multiple events per beam spill ob- 
served in the near-detector, the rate measured in the 
far-detector is approximately two events per day within 
the beam spill times, so the appropriate requirements for 
event selection are necessarily different. Specifically, the 
probability that two or more neutrinos produced in the 
same 10 fis beam spill will interact in the far-dctcctor is 
negligible. Therefore, if two events are reconstructed in 
the same spill, the coincidence is due either to a recon- 
struction failure or else one of the events has a nonbeam 
origin. Effects of multiple event reconstruction are miti- 
gated by requiring an event to be used in the analysis to 
contain at least 75% of the total deposited energy during 
the beam spill. 

The main background in the far-detector results from 
detector noise arising from the electronics and PMTs or 
from spontaneous light emission from the scintillator and 
wavelength-shifting fibers [13] ■ The noise from the elec- 
tronics and PMTs is removed by setting an energy thresh- 
old in the PMTs. The spontaneous light emission is re- 
moved by requiring accepted events to include at least 
nine strips or at least 10MIP deposited in the detector. 
Alternatively, events are also accepted if they include 
more than five strips and deposit more than 5MIP in 
the detector. 

Muons from cosmic rays are a potential source of back- 
ground events. Given the 0.2 Hz cosmic-ray muon rate 
at the far-detector, the number of cosmic-ray-induced- 
muons that may potentially coincide with beam spills is 
comparable to the number of beam-induced neutrino in- 
teractions observed. Most cosmic-ray-induced muons are 
well reconstructed and arc efficiently removed by the fidu- 
cial requirement. However, the reconstruction algorithms 
are optimized to handle recorded energy flow in the gen- 
eral direction of the beam. Problem cases can thus arise 
with very steep cosmic muons, which are removed by re- 
quiring the absolute value of the muon direction cosine 
in the longitudinal direction, \p z \/E, to be higher than 
0.4. In some cases the events are so steep that they are 
reconstructed only as showers and may be removed by 
using selection variables that describe the transverse and 
longitudinal shower profiles. The transverse variable is 
defined by calculating the root-mean-square (rms) value 
of the shower strip positions, whereas the longitudinal 
variable is defined as the ratio of active strips per plane 
to the total number of active planes in the event. Only 
those showers with a transverse rms value lower than 0.5 
and (strips/plane)/ (event planes) < 1 arc accepted for 



further analysis. Cosmic-ray muons that stop in the de- 
tector can mimic beam events if the end of the stopping 
muon track is interpreted as the vertex and the track 
is then reconstructed backwards. These events can be 
identified by performing a linear fit to the timing dis- 
tribution for strips on a track as a function of the strip 
longitudinal position. A fit resulting in a negative slope 
indicates that the event is a downward-going cosmic-ray 
muon and not a beam neutrino. The sample contamina- 
tion from cosmic-ray induced muons after these criteria 
are applied is estimated to be less than 0.1% [28j . 

Another potential background arises from data 
recorded while the Light Injection calibration system (LI) 
is flashing during normal data taking. The light injec- 
tion events are removed with 99.99% efficiency by using 
information from a PMT directly connected to the light 
injection system. By applying additional requirements 
based on concentrated detector activity, it is estimated 
that much less than one LI event is accepted in the entire 
data sample [28j . Furthermore, application of the LI re- 
jection criteria results in no measurable loss of efficiency 
for beam-neutrino interactions. 



IV. EVENT CLASSIFICATION 

After the selection criteria described in the previous 
section arc applied, the analysis proceeds by distinguish- 
ing charged-current events from neutral-current events. 
Distinct event classification procedures are employed for 
each sample, as described below. The reconstructed neu- 
trino energy spectra for both event classes are used in 
the fits described in Sec. I Villi and Sec. IIXI 

The goal of the event classification is to maximize 
the efficiency and purity of selected samples of neutral- 
current and charged-current events. Using Monte Carlo 
event samples, efficiency is defined as the number of true 
events of one type which are classified as that type, di- 
vided by the total number of true events of that type that 
pass the criteria described in Sec. Mil Purity is defined as 
the ratio of the number of true events of one type selected 
to the total number of events selected as that type. 

To avoid biases, the methods for identifying neutral- 
current candidate events and procedures employed in pre- 
dicting the far-detector spectrum, described in Sec. E] 
were developed and tested using only the near-detector 
data and Monte Carlo simulation. All analysis proce- 
dures were finalized prior to examining data in the far- 
detector. 

The neutral-current event classification employs sev- 
eral criteria based on reconstruction variables displaying 
large differences between neutral-current and charged- 
current events |29|. Charged-current events with short 
or no apparent tracks and poorly reconstructed events 
are the two main sources of background. The latter is 
mitigated by employing the various selections described 
in Sec. IIII1 The classification variables considered arc: 
event length, expressed as the difference between the first 
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and last active plane in the event; number of tracks re- 
constructed in the event; and track extension, defined as 
the difference between track length and shower length. 

A sample of candidate neutral-current events is ob- 
tained by applying specific requirements on the classifi- 
cation variables. Since neutral-current events are typi- 
cally shorter than charged-current events, events cross- 
ing fewer than 60 planes and for which no track is re- 
constructed are classified as neutral-current. Because 
neutral-current events are expected to have short or no 
reconstructed tracks, events crossing fewer than 60 planes 
that contain a track are classified as neutral-current if the 
track extends fewer than 5 planes beyond the shower. 
The values chosen maximize sensitivity for detection of 
sterile-neutrino admixture. Finally, events that are not 
classified as neutral-current-like are labeled as charged- 
current-like if they pass the classification procedures de- 
scribed in a previous MINOS publication 15 . These re- 
quirements arc applied to both near and far-detectors to 
obtain neutral-current and charged-current event sam- 
ples. 

The 6% of near-detector events and 3.5% of far- 
detector events that are not classified as either neutral- 
current or charged-current are not used in further anal- 
ysis. Evaluation of these removed samples using simu- 
lated data shows that approximately 75% of the events 
in each detector are charged-current interactions, and the 
reconstructed energies of those neutrinos distribute in ac- 
cord with the expectation based on the simulation. The 
remaining events are neutral-current interactions whose 
distributions in reconstructed energy are very similar in 
the two detectors; the latter distributions span the full 
visible energy spectrum, but with modest accentuation of 
the lower energy region GcV < ii^cco < 4 GcV. There- 
fore, the removal of these events does not introduce anal- 
ysis biases. 

Distributions for the event length and track extension 
classification variables for data of the near-detector are 
shown in Figs.[3ji and|3jD. The data are plotted together 
with the prediction of the MINOS Monte Carlo simu- 
lation, which adequately reproduces the shapes of the 
classification-variable distributions. Comparisons of dis- 
tributions in the far-detector for the same event classi- 
fication variables are shown in Figs. @Ji and |3Jd. Here, 
the Monte Carlo simulation uses oscillation parameter 
values obtained from the most recent MINOS charged- 
current measurement, |Ato3 2 | = 2.43 x 10 _3 eV 2 and 
sin 2 26» 23 = 1 [13. 

Efficiencies and purities for the neutral-current and 



b Candidate charged-current events are selected using a likelihood- 
based particle identification parameter constructed from three 
probability density functions. The functions represent distribu- 
tions for the variables (i) event length, (ii) fraction of the total 
event signal in the reconstructed track, and (iii) average signal 
per plane induced by the track. Further details are presented in 
Ref. 0. 
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FIG. 3: Comparisons of near-detector data with Monte Carlo 
predictions for distributions of the variables (a) event length 
and (b) track extension. The data quality requirements de- 
scribed in Sec. 1 1 1 1 1 are applied. Systematic uncertainties are 
displayed as shaded bands on the Monte Carlo expectation. 
Events selected as neutral-current-like are indicated by the 
arrows. 



charged-current event samples for both detectors are dis- 
played in Figs. [5] and [6] The classified neutral-current 
samples have nearly identical purities for the near and 
far-detectors. The far-detector purity curve' is computed 
for the case of null neutrino oscillations. The neutral- 
current sample efficiencies have identical trends but ex- 
hibit a constant relative offset over the full range of re- 
constructed event energies, E reco , due to the special near- 
detector selection criteria. The minima observed in sam- 
ple purities for both detectors near the peak of the fo- 
cussed neutrino spectrum, -E reC o ~ 4 GeV, reflect the 
abundant rate and consequently increased background 
from charged-current events. 

The resulting reconstructed energy spectra for the 
neutral-current and charged-current samples in the near- 
detector are shown in Figs. [7] and [8] respectively. The 
figures show that the selected neutral-current sample in- 
cludes a sizable background contribution from misiden- 
tificd charged-current events, in contrast to the se- 
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FIG. 4: Far-detector data versus predictions from the Monte 
Carlo (MC) simulation including — > v T oscillations, for 
distributions of the classification variables (a) event length 
and (b) track extension. The data quality requirements of 
Sec. IIIII are applied. The shaded bands show the systematic 
errors on the MC predictions. The arrows indicate events 
selected as neutral-current-like. 



lected charged-current sample which contains very little 
neutral-current background. For both samples, the data 
points are seen to fall within or near the 1 standard de- 
viation range of the systematic uncertainty of the Monte 
Carlo simulation. 
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FIG. 5: Selection efficiency and sample purity for Monte 
Carlo (MC) events selected as neutral-current candidates in 
the near and far-detectors, as a function of reconstructed 
event energy. Detector selection, fiducial volume and neutral- 
current/charged-current separation requirements have been 
applied. The shapes of both efficiency and purity distribu- 
tions are observed to be very similar for each of the two MI- 
NOS detectors. 
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FIG. 6: Distributions of selection efficiency (lower curves) and 
sample purity (upper curves) versus reconstructed event en- 
ergy, for Monte Carlo (MC) events selected as charged-current 
candidate events in the near-detector and in the far-detector. 



V. FAR-DETECTOR PREDICTION 



The predictions of the energy spectra of the neutral- 
current and charged-current samples at the far-detector 
are based on the observed near-detector data and make 
use of the expected relationship between the neutrino 
fluxes at the two sites. The process of making the pre- 
dictions is called "extrapolation" and may be viewed as 
making corrections to the simulation of interactions in 
the far-detector based on the energy spectrum measured 
in the near-detector. 

This analysis uses an extrapolation technique called 



the "far over near" (F/N) method [J|30|. This method 
makes the prediction of the far-detector spectrum by tak- 
ing the product of two quantities. The first quantity 
is the ratio of the expected number of events from the 
Monte Carlo simulation for each energy bin in the far- 
detector and near-detector spectra. The expected num- 
ber of events for each energy bin in the far-detector spec- 
tra depends on the composition of event types entering 
the samples and the corresponding oscillation probabili- 
ties for that energy. The second quantity is the number 



9 



40 



Near Detector Data 
Monte Carlo Expectation 
CC Background 




20 



(GeV) 



FIG. 7: Distribution of reconstructed visible energy for se- 
lected neutral-current events in the near-detector, for the data 
(solid points) versus the Monte Carlo prediction (open his- 
togram). The systematic errors (la) for the Monte Carlo are 
shown by the shaded band. Also shown is the Monte Carlo 
prediction for the background of misidentified charged-current 
events in the near-detector sample (hatched histogram). 
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FIG. 8: The reconstructed energy spectrum for selected 
charged-current events of the near-detector, for the data (solid 
points) versus the Monte Carlo (open histogram). Systematic 
errors at la for the Monte Carlo are indicated by the shaded 
band; the hatched histogram (lower left) shows the Monte 
Carlo expectation for the small background contribution by 
misidentified neutral-current events to this sample. 



of observed near-detector data events. The F/N method 
prediction is robust against distortions arising from dif- 
ferences between data and Monte Carlo simulation in the 
near-detector as these distortions are translated to the 
far-detector and do not affect the oscillation measure- 
ment 3. 

For example, the extrapolation for the neutral-current 
spectrum accounts for both the signal neutral-current 
events and the background charged-current events from 
each neutrino flavor. The extrapolation also addresses 
the ways in which the backgrounds change between the 



near and far-detectors. Specifically, for the case of the 
Vfi charged-current component of the neutral-current and 
charged-current samples, the F/N extrapolation predicts 
the number of events at the far-detector for the i-th bin 
of reconstructed energy to be 



^prediction ^ydata 



V 



(1) 



where is the number of selected events in the i-th 

reconstructed energy bin in the near- detector and Nf* 
is the number of events expected in that bin from the 
near-detector Monte Carlo simulation. The F^ c repre- 
sents the number of events expected from the far-detector 
Monte Carlo simulation in the i-th bin of reconstructed 
energy and j-th bin of true neutrino energy. In the equa- 
tion, Ej is the true neutrino energy and P Ull ^ Ux the prob- 
ability of muon-ncutrino transition to any other flavor. 

In particular, for the neutral-current spectrum, the 
extrapolation must take neutrino oscillations into ac- 
count to properly characterize the predominant back- 
ground arising from misidentified charged-current v^, 
and it must include the small spectral distortion re- 
sulting from misidentified charged-current v r and v e 
events. Thus, there are five separate classes of events 
that must be extrapolated to the far-detector: (i) genuine 
neutral-current interactions, (ii) charged-current in- 
teractions, (iii) v T charged-current interactions, (iv) pos- 
sible v e charged-current interactions originating from 
Vfi oscillations, and (v) charged-current v e interactions 
initiated by the intrinsic v e beam component. The muon 
neutrinos in the simulation include oscillations and are 
integrated in bins of reconstructed energy to account for 
the changing background. Oscillations of the intrinsic 
beam v e into are not taken into account as those v e 
comprise only 1.3% of the neutrinos in the beam and 
mixing angle for such oscillations is so small as to make 
the contributions from that oscillation mode negligible. 

The five classes of simulated events are used to con- 
struct individual two-dimensional histograms of true neu- 
trino energy versus reconstructed energy. In each of these 
histograms, all the events in an individual bin of true en- 
ergy are multiplied by the same survival or transition 
probability. After applying oscillations the simulated 
events are converted to a reconstructed energy spectrum 
by integrating across all the true energy bins for each bin 
of reconstructed energy, producing the sum in Eq. (JTJ) . 
The reconstructed energy spectra for each separate data 
set, signal and backgrounds, are added together into one 
spectrum. 



VI. SYSTEMATIC UNCERTAINTIES 

The principal sources of systematic uncertainties in 
these analyses are: (i) absolute scale of the hadronic en- 
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ergy; (ii) relative calibration of the hadronic energy be- 
tween the two detectors; (iii) relative normalization be- 
tween the two detectors; (iv) charged- current background 
in selected neutral-current events; and (v) uncertainties 
due to the near-detector selection requirements in the 
near-detector event counts. Monte Carlo studies have 
been performed where, for each single uncertainty, the 
Monte Carlo spectrum is varied by ±1 standard devia- 
tion independently, in order to estimate the effect of each 
on the extrapolated spectrum. Beam and cross section 
uncertainties that are common to the two detectors ef- 
fectively cancel when using the F/N extrapolation. 

The absolute hadronic energy scale has an uncertainty 
of 12%. This value is a combination of the uncer- 
tainty in the hadronization model and intranuclear ef- 
fects (10%) and uncertainty of the detector response to 
single hadrons (6%) . An uncertainty of 3% on the relative 
energy scale between the two detectors was determined 
from cross-calibration studies using stopping muons 
The relative normalization between the two detectors has 
an uncertainty of 4%. This is a combination of the uncer- 
tainties due to fiducial mass, live time, and reconstruc- 
tion differences between the two detectors. 

To evaluate the uncertainties due to the near-detector 
selection, the requirement that the total number of recon- 
structed strips in an event is at least four was shifted by 
±1 strip. The effects on the reconstructed energy spec- 
trum were determined for each shift. The uncertainty has 
been estimated to be 15.2% for E rcco < 0.5 GeV, 2.9% 
for 0.5 < E ICCO < 1 GeV, 0.4% for 1 < E ICCO < 1.5 GcV 
and is negligible for higher visible energies. 

The uncertainty in the number of charged-current 
background events is determined using near-detector 
data taken in several different beam configurations, 
namely, (i) horns-off, in which there was no current in the 
magnetic horns; (ii) medium energy, in which the target 
is moved upstream from the first horn by 100 cm; and 
(iii) high energy, in which the target is moved upstream 
from the first horn by 250 cm. For each of these beam 
configurations, described in further detail in Ref. [H, the 
charged-current background component has a consider- 
ably different spectrum from the one obtained in the low- 
energy (LE) configuration. The charged-current back- 
ground component in the neutral-current spectrum can 
thus be determined by using the observed differences in 
energy spectrum between the low-energy beam configura- 
tion and each of the other beam configurations along with 
information from Monte Carlo simulation of each config- 
uration. In the low-energy configuration, the number 
of selected near-detector neutral-current events JV LE can 
be written as the sum of true neutral-current events and 
background charged-current events in that beam config- 
uration: 



of selected neutral-current events may be written as 



A LE = N% + A C LE 



(2) 



N Mt = r Mt ■ N hE -4- r Alt • /V LE 
nc 1 v nc ' cc cc i 



where = N^/N^f and r, 
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(3) 



are 



determined from the Monte Carlo simulation. Equa- 
tions ([2]) and <j3j> can be solved to yield the solutions: 
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(4) 



The final estimate on the iV^. E background number re- 
sults from the weighted average of the three different so- 
lutions of Eq. (|4]) when the LE beam data is compared 
with each of the other beam configurations. The un- 
certainty on the N^r background number is 12% ± 2%. 
Therefore, the uncertainty in the number of charged- 
current background events is conservatively taken to be 
15% at all energies at the near and far-detectors. 



VII. SEARCH FOR ACTIVE NEUTRINO 
DISAPPEARANCE 

The data collected in the near and far-detectors have 
been classified as either neutral-current or charged- 
current events using the selections described in the pre- 
vious sections. A total of 388 data events are selected as 
neutral-current in the far-detector. The measured and 
predicted -Ereco spectra at the far-detector are shown in 
Fig. The prediction assumes that oscillations occur 
only among the three active flavors and uses the values 
of | Attt.§ 2 I and 6*23 previously measured by MINOS [Io| . 

Although the present analysis is not capable of isolat- 
ing an electron neutrino appearance signal, it must take 
— > v e oscillations into account because the classifica- 
tion criteria of this analysis include v e charged-current 
interactions in the neutral-current enriched sample with 
nearly 100% efficiency. This accounting is done by com- 
paring the observed neutral-current spectrum to two pre- 
dictions, one that assumes null v e appearance, and an- 
other that assumes an upper limit for the v e appear- 
ance rate in the far-detector calculated with the normal 
neutrino-mass hierarchy, 6*13 = 12° and 6 = 3ir/2. The 
choice of 6*13 corresponds to the 90% confidence level 
upper limit established by the CHOOZ reactor experi- 
ment [U for the |Am§ 2 | value measured by MINOS [Toj ] ; 
the choice of S maximizes the v e appearance for the cho- 
sen value of 6*13. As seen in Fig. [SI the observed spectrum 



matches the prediction based on oscillations among the 
three active flavors very well over the full range of allowed 
values of #13. 

The agreement between the observed and predicted 
neutral-current spectra is quantified using a statistic, R: 



B 



Apata — B C C 

Snc 



(5) 



For an alternate beam configuration "Alt" , the number 



where, within a given energy range, Apata is the observed 
event count, i?cc is the extrapolated charged-current 
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4 — Far Detector Data 




E reco (GeV) 

FIG. 9: The reconstructed energy spectrum ol neutral-current 
selected events at the far-detector (points with statistical un- 
certainties). The Monte Carlo prediction assuming standard 
three-flavor oscillations is also shown, both with (dashed line) 
and without (solid line) i/ e appearance at the CHOOZ limit. 
The shaded region indicates the 1 standard deviation system- 
atic uncertainty on the prediction. The hatched region shows 
the Monte Carlo prediction for the background of misidenti- 
fied charged-current events in this sample. 

background from all flavors, and Snc is the extrapolated 
number of neutral-current interactions [l"7| . The values 
of Snc and contributions to Bqc are shown in Table |U 

The disappearance of occurs mainly for true neu- 
trino energies < 6 GeV [10| . While the true energy of 
a neutrino interacting through the neutral-current pro- 
cess cannot be measured, the data can be separated into 
two samples whose reconstructed energy roughly discrim- 
inates between neutrinos with true energies greater than 
and less than 6 GeV. According to the Monte Carlo 
simulation, events with E ICCO < 3 GeV have a me- 
dian true neutrino energy of 3.1 GeV while events with 
3 GeV < E TfiCO < 120 GeV have a median true neutrino 
energy of 7.9 GeV. The values of R calculated for these 
ranges in E rcco are shown in TablelU For all data samples, 
R differs from unity by less than 1 standard deviation. 
The effect of v e appearance on the value of R is treated 
as an uncertainty in this analysis and is indicated by the 
last uncertainty in Table [J] Because j/ e -charged-current 
events are almost always identified as neutral-current, the 
effect of v e appearance at the far-detector is to decrease 
the value of i?, since the predicted background will be 
larger than for the case of null v e appearance. 

The measured values of R for each energy range in- 
dicate that neutrino oscillations among the active fla- 
vors describes the observed data quite well. Over the 
full energy range, — 120 GeV, a value of R = 1.04 ± 
0.08(stat.) ± 0.07(syst.) — 0.10(^ e ) is measured, corre- 
sponding to a depletion of the total neutral-current event 



rate assuming null (maximally-allowed) v e appearance of 
less than 8% (18%) at 90% confidence level. The follow- 
in g sections explore the extent to which the data allo w 



-E rcco (GeV) iV D ata <Snc -^cc 


n cc 


B v e 
cc 


- 


3 


141 


125.1 13.3 


1.4 


2.3 (12.4) 


3- 


120 


247 


130.4 84.0 


4.9 


16.0 (32.8) 


- 


3 


R = 


0.99 ±0.09 ±0.07 


-0.08 




3 - 


120 


R = 


1.09 ±0.12 ±0.10 


-0.13 




- 


120 


R = 


1.04 ±0.08 ±0.07 


-0.10 





TABLE I: Values of R, iVoata, SWc, an d the contributions to 
Sec for various reconstructed energy ranges. The numbers 
in parentheses are calculated including u e appearance at the 
CHOOZ limit, as discussed in the text. The first uncertainty 
in the value of R shown is the statistical uncertainty, the sec- 
ond is systematic uncertainty, and the third is due to possible 
v e appearance. 

oscillations between the active flavors and one sterile neu- 
trino or oscillations in combination with neutrino decay. 



VIII. OSCILLATIONS INCLUDING A STERILE 
NEUTRINO FLAVOR 

Mixing of the three active neutrino flavors with one 
sterile neutrino requires the addition of one mass eigen- 
state. The probability for mixing between any two flavors 
is described by the neutrino mixing matrix [32[, U, which 
defines the rotation from the mass basis into the flavor 
basis. Thus, the mixing matrix needs to be expanded by 
one row and one column to accommodate the additional 
mass and flavor states. The expanded 4x4 mixing ma- 
trix contains six mixing angles and six phases, with three 
of the phases being Majorana phases that are not rele- 
vant to an oscillation experiment [33j . The matrix can be 
written as a product of six independent rotation matri- 
ces about the Euler axes Rij , where i j refers to the plane 
in which a particular rotation takes place. The ordering 
of the rotation matrices is arbitrary, reflecting that the 
mixing matrix can be parameterized in many ways. The 
ordering described below was chosen to make the analysis 
more straightforward. 

MINOS is designed to precisely measure |Am§ 2 | but 
has no sensitivity to Am^. Consequently, the mass 
states 1 and 2 are treated as degenerate in this analy- 
sis. When two mass states are degenerate the rotation in 
the ij-plane is unphysical and the corresponding mixing 
angle, 6ij, vanishes from the oscillation probabilities [l2T |. 
Additionally, the matrix should be of the form such that 
the U e 3 component of the mixing matrix goes to zero 
when 6*i3 = to distinguish an electron component in 
the third mass eigenstate from effects of sterile neutrinos. 
For those reasons, the general form of the mixing matrix 
used by the current analysis is written as 
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U — i?34(^34)^24(^24, ^2)^14(^14)^23(^23)^13(^13, <5l)-Rl2(^12, <53) 
= i?34(^34)^24(^24, ^2)^14(^14)^23(^23)^13(^13, <5l), 

I 



(6) 



where the 5k are the three CP-violating Dirac phases and 
the last equality reflects the assumption of degeneracy in 



J 



mass states 1 and 2. Thus, the mixing matrix can be 
written as, 



U 



U T2 
U s2 



ci4,si3e~ lSl 
-s 14 s 13 e- lSl S24e~' lS2 +Ci 3 s 23 C24 

-Si4C24S34Sl3e~ tl51 - Ci 3 S2 3 S34S24e 4 ' 52 + C13C23C34 
-Sl4C24C34Sl3e~ i<51 - Ci 3 S 2 3C 3i S 2 4e l52 ~ C13C23S34 



l& 2 



su 
c 14 s 2 4e 
C14C24S34 
C14C24C34 



(7) 



r 



Here Cy = cos 



?ij and Sij = sin ( 



ij , and only the elements 
of the matrix that appear in the oscillation probabilities 
given below have been expressed in terms of the mixing 
angles and phases. 



In this equation U^j is the element of the mixing matrix 
describing the coupling between the muon flavor state 
and mass state j. By squaring the amplitudes and using 
the unitarity of U, one obtains the oscillation probabili- 
ties for the different channels, 



A. Oscillation Probabilities 



The oscillation probabilities are derived following the 
example in Ref. [341 ] . The evolution of a neutrino with 
initial flavor state a is given by 



i^) = E^ e " i; " lj2L/(2£) i^'°) 
3 



(8) 



where the sum is over the mass states, E is the neutrino 
energy, L is the distance traveled by the neutrino and mj 
is the mass of state j. As the NuMI beam is almost en- 
tirely or V/j,, only the amplitudes and probabilities for 
Vp — > u x , where x represents e,/x, r, or s, are described. 
The amplitude for — > v x is given by 



-im?L/(2E) 



(9) 



P(u a -y up) = 6 a ?-4j2 n (Kj U lH u cAU}i) sin 2 A Tl 
+2 AKj UpiUdVti) sin 2A ji , (10) 



where Aji = (m| — m 2 )L/(4:E) and 1ZQ, X() designate 
the real and imaginary parts of the amplitude products. 
Because A21 <C A31 the first and second mass states are 
treated as degenerate and the factors of SU1A21 can be 
set to 0. This degeneracy also implies that A 4 2 = A41 
and A32 = A31 . Equation (fTU)) can be expanded for the 
different oscillation scenarios, 



Pu^u, =1-4 |C/ M3 | 2 (1 - \U^\ 2 - IM 2 )sin 2 A31 + |t/ M 4| 2 |t^3| 2 sin 2 A 43 + \U^\ 2 {1 - {U^ 2 - \U^\ 2 ) sin 2 A41 



ATll \U^\ 2 \U a3 \ 2 sin 2 A31 + \U^\ Z \U ai \ 2 sin 2 A41 + U*^U aA U&U^{saf A 31 - sin 2 A 43 + sin 2 A 



+2l|L/; 4 L/ Q 4C/ AI 3f/: 3 (sin2A3i - sin2A 41 + sin2A 43 )|, 



(11) 



r 



where a = e, r, or s and the orthogonality constraint mcnts corresponding to the first and second mass states. 
J2i UaiUfo — has been used to eliminate the matrix ele- 
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FIG. 10: Schematic representation of the possible mass spec- 
tra for models including four mass eigenstates. In the left 
panel, the first and fourth eigenstates are degenerate, the cen- 
ter panel has the fourth mass eigenstate much heavier than 
the third, and the right panel has the third and fourth eigen- 
states as degenerate. Only the scenarios illustrated in the left 
and center panels are considered in this analysis. 



As can be seen from Eqs. (JT]) and (|TT|) . mixing between 
three active and one sterile neutrino in the most general 
case involves ten parameters which are mostly unknown: 
five mixing angles, three mass-squared splittings, and two 
CP-violating phases. The following simplifying assump- 
tions are made to reduce the number of possible parame- 
ters. First, 0i3 is eliminated as a free parameter by only 
considering two values, $13 = 0° or #13 = 12°, where the 
latter is the CHOOZ limit at the MINOS measured value 
of |Am| 2 |. When 0i 3 is set to the CHOOZ limit, S 1 is 
taken to be 37r/2 as that is the value of the CP- violating 
phase that maximizes the v e appearance probability. Sec- 
ond, the additional CP- violating phase is removed by set- 
ting 82 = as MINOS has no sensitivity to this phase. 
Finally, nonzero values of #14 do not measurably change 
the oscillation probabilities, so that angle is set to 0° as 
well. 



In addition to the parameters represented in the mix- 
ing matrix and oscillation probability equations, there is 
one more free parameter, namely the mass of U4 rela- 
tive to the other mass states. To limit the number of 
free parameters in the models, this analysis only exam- 
ines possibilities for the relative mass of 1/4 that allow the 
oscillation probabilities to depend only on lAm^j. The 
possible hierarchies are illustrated in Fig. [TU] and only 
the normal mass hierarchy is shown. An inverted mass 
hierarchy would change the relative locations of the pair 
V\ and v% to ^3. The fourth mass state could either be 
degenerate with the first mass state, as seen in the left 
panel of the figure, much more massive than the third 
mass state, as seen in the middle panel, or degenerate 
with the third mass state, as shown in the right panel. In 
the scenario with the third and fourth eigenstates being 
nearly degenerate there would be no discernible mixing 



between the active and sterile components. This conse- 
quence comes from the SNO results [f| which indicate 
that any coupling between v s and the active neutrinos 
occurs in the third and fourth mass eigenstates. There- 
fore, in this case the mass difference in which oscillations 
between active and sterile neutrinos could arise is too 
small to be observed in terrestrial baselines and it is not 
explored further in this analysis. 

An alteration of the oscillation probability of Eq. (fTTj) 
is predicted for neutrinos which traverse dense matter. 
By virtue of their neutral-current coherent forward scat- 
tering with nucleons, active neutrinos acquire an effective 
matter potential whereas sterile neutrinos do not. The 
effective potentials arc identical for the and v T fla- 



vors so that matter effects vanish for v u 



mixing. 



For — > v s oscillations, however, the overall nonzero 
matter potential yields MSW-type modifications to the 
mixing angle and oscillation length such that, in normal 
(inverted) hierarchy scenarios, the oscillation probab ility 
will be suppressed (enhanced) relative to — > v r [351 ] . 
Observations of neutrinos with energies above 12 GeV 
that have traveled through several thousand kilometers 
of dense material, such as those described in Ref. [||, 
would be sensitive to this matter effect. By contrast, the 
beam-induced neutrinos monitored at the MINOS far site 
originate from a E„ spectrum highly peaked between 2 
and 6 GeV and travel 735 km through the Earth's crust. 
For the effective mixing angles and oscillation lengths of 
neutrinos in the NuMI beam, the — > v s matter effect 
only gives rise to subpercent distortions of the neutrino 
oscillation probability that are mostly confined to neu- 
trino energies below 2 GeV. Distortions of this magnitude 
are of no consequence for the studies reported here, and 
so matter effects involving sterile neutrinos have been 
neglected. 



1. Active-sterile mixing when 7714 = mi 

In the m.4 = mi case, the first and fourth mass eigen- 
states are assumed to be degenerate. Because the first 
and second eigenstates are also treated as degenerate, the 
second and fourth states are degenerate as well. These 
degeneracies allow one to set 9u = 824 = 0° in Eq. (0), 
which reduces the number of parameters in the model to 
four. There are no v e or components in the fourth 
mass eigenstate, however there is a v s component in the 
third mass eigenstate. Using these simplifications, the 
oscillation probabilities become 



P 



_>„„ = 1-4|IV 3 | 2 l-lCU 2 )sin 2 A 31 , 



P v ^ Va = 4|C/ M3 | 2 |C/ Q3 | 2 sin 2 A31 



(12) 



This model is equivalent to the 
presented by MINOS in Ref. 



le p 



henomenological model 
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2. Active- sterile mixing when 7774 3> m-j 

In the TO4 3> "^3 case, the fourth mass eigenstate is as- 
sumed to be much larger than the third; consequently the 
values of sin 2 A 41 and sin 2 A 43 average to \ . Addition- 
ally, sin2A4i and sin2A 4 3 average to 0. In this model 
Am| 3 is assumed to be C(eV 2 ) such that the regime 
of rapid oscillations and thus the averages mentioned are 
valid at the far-detector site, while ensuring no detectable 
depletion of occurs at the near-detector. Such models 
have recently received attention in the literature [HI, [l3| . 
Using the above simplifications reduces the number of 
parameters in this model by two, and allows Eq. (fTTj) to 
be written as 



i - ^ ic/ m3 | 2 ( i - IM* - IM* ) sin " A 3i 



+ J ^f-(i-|£U') 

A%\ [ \U^\ 2 \U a3 \ 2 + U^U^sKs ) sin 2 A31 



\Ufi4 \ '\Vai\ 



(13) 



where the second term of Eq. (fTTj) does not appear be- 
cause of the assumptions that #14 = 0° and 62 = 0. 



B. Fitting active-sterile oscillations to the data 

The data are compared to Monte Carlo predictions 
based on the probabilities in Eqs. (Tl"2"|) and ([To]) using 
the x 2 statistic appropriate to small sample sizes, 



N 



Oi In — 



.7=1 3 



(14) 



Here ej is the expected number of events, assuming os- 
cillations among four flavors, in bin i of the energy spec- 
trum, and Oi is the observed number of events in that 
bin. The second sum is the contribution to x 2 from the 
parameters describing the systematic uncertainties; the 
nuisance parameter Ej is the shift from the nominal fit 
value for the j-th source of systematic uncertainty and 
(Tj is the uncertainty associated with that source. Both 
the neutral-current spectrum shown in Fig. [9] and the 
charged-current spectrum shown in Fig. [TT] arc used to 
obtain the oscillation parameters that best fit the data. 
The neutral-current spectrum provides information on 
the mixing angles for mixing between active and sterile 
neutrinos while the charged-current spectrum provides 
constraints on the mixing angle #23 and the mass splitting 
Am 2 .^. The charged-current-like spectrum of Fig. [TTlis 
statistically consistent with that presented in Ref. [10| 
given the different fiducial volumes and event separation 
procedures used in the two analyses. The five system- 
atic uncertainties described in Sec. IVII arc included as 
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FIG. 11: The reconstructed energy spectrum of charged- 
current selected events at the far-detector (points with sta- 
tistical uncertainties). The Monte Carlo prediction assuming 
standard three-flavor oscillations (solid line) is shown with 
the 1 standard deviation systematic uncertainty on the pre- 
diction indicated by the shaded region. The prediction for 
the small background of misidentified neutral-current events 
in this sample (hatched region) is also shown. 



nuisance parameters in the far-detector fits. By fitting 
for the systematic parameters simultaneously using both 
near and far-detector data, the effect of the uncertainties 
is substantially reduced due to significant cancellations 
of uncertainties between the two detectors. 

The best-fit values for the mixing angles in the two 
models as well as the % 2 for each are shown in Table [TTJ 

A list of the systematic effects for the fit parameters 
is presented in Table IIIII The latter table shows that for 
each mixing angle evaluated in the four-neutrino models 
considered, the uncertainties introduced by the five most 
significant sources of systematic error are relatively small 
compared to the total uncertainty ranges obtained from 
the fits, as summarized in Table [TT1 

The best-fit values obtained for each model for the 
mass splitting, |Am 31 |, agree well with the result found 
in Ref. [Toj | . The one-dimensional projections of the A% 2 



Model 


013 


x 2 /d.o.f. 


623 


#24 


634 


/. 


7714 = 77ll 



12 


47.5/39 
46.2/39 


45.0l 9 8 ;° 
47 1 + 8 - 8 




i+ 28 - 7 

23.0l24.i 


0.51 
0.55 


7774 2> 7773 



12 


47.5/38 
46.2/38 


45.0l|;° 

47 1 +8.8 


o+ 7 - 2 

u - u -o.o 

n o+ 7 - 2 

u - u -o.o 


1+ 28 ' 7 
U -- L — 0.1 

23.0l 22 ;i 


0.52 
0.55 



TABLE II: Best-fit points and uncertainty ranges obtained 
for the active-sterile oscillation models. Results are shown 
with and without i/ e appearance at the CHOOZ limit. All 
angles are given in degrees. The quantity f s is defined as the 
fraction of disappearing that could transition to v a and is 
given at the 90% C.L. in this table. The values of f a in the 
7774 2> 7713 model are evaluated for E v = 1.4 GeV. 
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Model 


Parameter 


Shift due to Systematic Uncertainty 






Absolute Enad. 


Relative Enad. 


Normalization 


CC Background 


ND Selection 


Total 




A a 
A6>23 


U.O 


U.O 


U.O 


A 1 
U.l 


U.l 


U. r 


7714 = TUl 
















iA(734 


3.6 


9.9 


12.6 


9.9 


9.9 


21.6 




At?23 


U.z 


U.O 


U. 1 


n o 
u.z 


U.z 


u. f 


m 4 > m 3 


Aft,,, 


1.5 


2.1 


5.1 


0.3 


0.3 


5.7 




A6> 3 4 


4.5 


9.9 


6.3 


9.9 


9.9 


18.8 


Oscillations 


Aq (GeV/km) 


2.54 x 10" 4 


0.70 x 10" 4 


6.25 x 10" 4 


1.23 x 10^ 4 


1.15 x 10" 4 


6.99 x 10" 4 


with decay 


A0 


2.6 


3.7 


0.9 


4.0 


3.9 


7.2 



TABLE III: Summary of mixing-angle deviations introduced by the major systematic uncertainties from best-fit results in 
which systematic shifts have been neglected. Angular deviations, shown in degrees, are displayed for each mixing angle fitted, 
for each of the neutrino models analyzed in this work. 
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FIG. 12: Projections of Ax 2 as a function of the mixing angles 
for the 77i4 = TTii model. The solid line contours are obtained 
with null v e appearance, whereas the dashed line contours 
include v e appearance at the CHOOZ limit. The ranges of 
values allowed at 68% and 90% confidence levels lie within 
contours below the horizontal dashed lines. 



between the best-fit point and the remaining points in the 
space are shown in Fig. [T^] for the 7714 = 7771 model. The 
two-dimensional 90% confidence level contours for that 
model arc shown in Fig. 1131 The projections and contours 
for the 7774 3> m.3 model are shown in Figs. [141 and [T5l 

As seen in these figures, #34 < 38° (56°) at 90% con- 
fidence level for the 7774 = m\ model. The number in 
parentheses represents the 90% confidence level limit ob- 
tained when maximal v e appearance is allowed. For the 
m 4 > m 3 model, 6» 24 < 10° (11°) and 6> 34 < 38° (56°) at 
the 90% confidence level. These limits indicate that any 
coupling between the active neutrinos and a sterile neu- 
trino is submaximal. Furthermore, the \ 2 values indicate 
that these four-flavor models fit the data no better than 
oscillations among only the active neutrinos. 

A straightforward method to quantify the coupling be- 
tween the active and sterile neutrinos is to determine the 
fraction of disappearing Un that transition to v s . That 
fraction is expressed as 



fs 



p v 



(15) 
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13: Contours representing 90% confidence level for the 
7774 = mi model. The solid line and best-fit point (solid sym- 
bol) are obtained assuming null v e appearance, whereas the 
dashed line and corresponding best-fit point (open symbol) 
are obtained with v e appearance set at the CHOOZ limit. 



For the m 4 = mi model, the disappearance fraction f s 
is energy independent, as can be seen upon inserting the 
expressions from Eq. (fT2|) into Eq. (|T5|) . The 90% confi- 
dence level limit for f s is determined by selecting a large 
number of test values of 623 and #34 from Gaussian distri- 
butions with mean and a given in Table [XT] The value of 
f s that is larger than 90% of the test cases represents the 
limit. The value corresponding to the 90% confidence 
level is f s < 0.51 (0.55), with the value in parentheses 
indicating the value obtained for maximally- allowed v e 
appearance in the beam. This new limit on the value of 
f 3 represents a reduction of 33% compared to the previ- 
ous MINOS result without v e appearance (l7l |. 

For the 7774 3> model, inserting the expressions 
from Eq. (fT3|) into (fT5|) shows that f s is energy depen- 
dent because of the constant terms in Eq. (|13|1 . For this 
reason the 90% confidence level value of f s for this model 
is presented at E v = 1.4 GcV, the energy where the 
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FIG. 14: Projections of Ax 2 as a function of the mixing angles for the 7714 3> mz model. The solid line is obtained for the 
case of null v £ appearance whereas the dashed line represents solutions with i/ e appearance at the CHOOZ limit. The ranges 
of values allowed at 68% and 90% confidence levels lie within contours below the horizontal dashed lines. 
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FIG. 15: Contours representing 90% confidence level for the 
7Ti4 2> 7T13 model. The solid line and best-fit point (solid sym- 
bol) are obtained for the case of null v a appearance, whereas 
the dashed line and corresponding best-fit point (open sym- 
bol) is obtained with appearance included with #13 at the 
CHOOZ limit. 



disappearance probability is a maximum. The determi- 
nation of the limit follows the procedure described above, 
but with the addition of selecting a value of #24 for each 
test case as well. At 90% confidence level f s < 0.52 (0.55) 
for E v = 1.4 GeV in this model. Thus, in either model, 
approximately 50% of the disappearing can convert to 
v s at 90% confidence level as long as the amount of v e ap- 
pearance is less than the limit presented by the CHOOZ 
collaboration. 



IX. OSCILLATIONS WITH DECAY 

It was noted more than a decade ago that neutrino 
decay, as an alternative or companion process to neu- 
trino oscillations, offers some capability for reproducing 
neutrino disappearance trends [18j |. The model investi- 
gated here [36[ includes neutrino oscillations occurring in 
parallel with neutrino decay. Normal neutrino-mass or- 
dering is assumed, and the mass eigenstates V\, V2 are 
approximately degenerate, so that 3> vn-i « mi. The 
heaviest neutrino-mass state 1/3 is allowed to decay into 
an invisible final state. With these assumptions, and ne- 
glecting the small contributions from v e mixing, only the 
two neutrino flavor states v ^ and v T , and the correspond- 
ing mass states and ^3, are considered. The evolution 
of the neutrino flavor states is given by [36[ : 



. dv 

dx 



4E 

. m 3 

■1 

4r,£ 



-cos 29 sin26> 

sin 29 cos 29 

2sin 2 6> sin 29 

sin26> 2cos 2 ( 



0, (16) 



where T3 is the lifetime of the ^3 mass state and 9 is the 
mixing angle governing oscillations between and v T . 
Solving Eq. (fTBj) one obtains probabilities for survival 
or decay: 



Pu 



decay 



cos 4 9 + sin 4 9e 



2 cos 2 9 sin 2 9e 2T i E cos 
1 - e^^) sin 2 



™3 L ( Aml,X 



'32 J 

2E 



(17) 
(18) 



The limits T3 — > 00 and Am| 2 - ► correspond to a pure 
oscillations or a pure decay scenario, respectively. 

In a conventional neutrino oscillations scenario, the ra- 
tio of the predicted charged-current spectrum in the far- 
detector with the null-oscillation expectation displays the 
characteristic "dip" at the assumed Am| 2 value that is 
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Model 


X 2 /D-O.F. 


a (GeV/km) 


e 


Osc. with Decay 


47.5/39 


O.OOlH x lO" 3 


45.0™ 


Pure Decay 


76.4/40 


4.6±i;£ x 10" 3 


50.9t 3 ?i 7 



TABLE IV: Best-fit points and uncertainty ranges obtained 
for the relevant parameters of the oscillation with decay 
model. The result obtained for the pure decay scenario, 
Am§ 2 — > 0, is also presented. Angles are shown in degrees. 
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FIG. 16: The best-fit point and 90% C.L. contour for the two 
parameters of the neutrino oscillations-with-decay model, the 
neutrino mass-lifetime ratio a and the mixing angle 0. 



absent in the equivalent ratio computed for pure neutrino 
decay. Previously published results by MINOS using the 
charged-current far-detector spectrum support conven- 
tional oscillations and disfavor a scenario of pure neutrino 
decay at 3.7 standard deviations [Tol] . In the present 
analysis, both the neutral-current and charged-current 
far-detector spectra shown in Figs. I§1 and ITT1 arc included 
in the fit. Consequently, additional sensitivity is gained 
with respect to previous analyses, since any neutrino de- 
cay into a noninteracting final state would also deplete 
the neutral-current spectrum according to Eq. (|18j) . Ad- 
ditionally, the analysis is extended to the more general 
scenario combining oscillations and the decay model de- 
scribed above. 

The best-fit values extracted for 9 and the parameter 
a = wis/t^ using this model are summarized in Table iPvH 
Figure [TBI shows the two-dimensional 90% confidence in- 
terval obtained by the fit. The results are consistent with 
maximal mixing (9 = 45°) and with no neutrino decay 
(a- 
Rcf. 

Figure [T7] shows the one-dimensional A% 2 projections 
for a and 9, with other parameters marginalized. The 
90% confidence level limit found for the neutrino decay 
lifetime is r 3 /m 3 > 2.1 x 10~ 12 s/cV. 

A Ax 2 of 28.9 is obtained for the pure decay scenario. 
Thus, a pure neutrino decay model with null oscillations, 
as considered in Ref. [l(|, is disfavored at the level of 



The best-fit value for |Am| 2 | is consistent with 
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FIG. 17: Projections in for fit solutions for a and 

mixing angle for the oscillations-with-decay model. Parame- 
ter ranges allowed at 68% and 90% confidence levels lie below 
the corresponding dashed horizontal lines. 



5.4 standard deviations. 



X. SUMMARY 

Searches for depletion or distortion in rate and visi- 
ble energy spectra of neutral-current events recorded in 
the MINOS far-detector have been carried out for the 
purpose of detecting or constraining processes involving 
active-sterile neutrino mixing, as well as further restrict- 
ing models including neutrino decay. The data exposure 
analyzed corresponds to 3.18 x 10 20 protons on target col- 
lected during the period of May 2005 to July 2007. 

A total number of 388 neutral-current events were 
observed in the far-detector, whereas the expecta- 
tion from standard three-flavor neutrino models is 
377±19.4(stat.)±18.5(syst.) events. The value for the 
statistic R that gauges the agreement between the data 
and the expectation based on oscillations among the 
three active flavors is R = 1.04±0.08(stat.)±0.07(syst.)- 
0.10(V e ), which is consistent with no depletion of the ac- 
tive neutrino flux. 

Joint fits to the observed neutral-current and charged- 
current energy spectra, assuming two different neutrino 
oscillation models that include an additional sterile neu- 
trino flavor, yield the following 90% confidence level 
limits on the oscillation parameters: #34 < 38° (56°) 
for the model with 7714 = m\\ #24 < 10° (11°) and 
6*34 < 38° (56°) for the model with m 4 > m 3 . The val- 
ues in parentheses represent the results for maximally 
allowed v e appearance. These limits for the mixing an- 
gles between active and sterile neutrinos show that mix- 
ing between the active flavors dominates oscillations. In 
fact, the fraction of active neutrinos that oscillate into 
a sterile species is constrained to be f s < 0.51 (0.55) at 
the 90% confidence level for the m,^ = vn\ model and 
f s < 0.52 (0.55) for E v = 1.4 GeV in the case of the 
model with TO4 3> ms- 

Similar fits, assuming a two-flavor neutrino model in 
which oscillations may occur in parallel with decay into 
a sterile species, yield the best-fit values for the mass- 
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lifetime ratio a = O.OOIqq x 10~ 3 GcV/km and mixing 
angle 9 = 45.0°ig° g g 3 . From these results, we extract a 
90% confidence level limit on the neutrino decay lifetime, 
> 2.1 x 10~ 12 s/cV. The pure decay scenario 
(Am| 2 — >0) is disfavored by 5.4 standard deviations in 
our data as an alternative explanation to neutrino oscil- 
lations. 
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